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(57) Abstract: Some embodiments of the invention provide a location -determination system that includes several transmitters and 
at least one receiver. Each transmitter transmits a signal that includes a unique periodically-repeating component, and the receiver 
receives a reference signal. Based on the received reference signal, the location -de termination system identifies an estimated location 
of the receiver as follows. For each transmitter in a set of transmitters, the system computes a phase offset between the received 
reference signal and a replica of the transmitter's periodically-repeating component. The system also identifies an approximate 
location of the receiver and an approximate receive time for the received signal. The system then uses the identified approximate 
location and time, and the computed phase offsets, to compute pseudoranges for the set of transmitters. Finally, the system identifies 
the estimated location of the receiver by using the computed pseudoranges. 
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Location-Determination Method and Apparatus 

FIELD OF THE INVENTION 

The invention is related to location-determination method and apparatus. 

BACKGROUND OF THE INVENTION 

In many applications, it is necessary to estimate the location of objects in their 
environment. To date, numerous location-determination systems have been proposed 
for this task. One such system is the global positioning system (GPS). This system 
includes a number of satellites that orbit the Earth. It also includes GPS receivers, 
monitoring stations, and differential GPS receivers on Earth. 

GPS satellites transmit signals from which GPS receivers can estimate their 
locations on Earth. A GPS satellite signal typically includes a composition of: 
»(1) carrier signals, (2) pseudorandom noise (PRN) codes, and (3) navigation data 
GPS satellites transmit at two carrier frequencies. The first carrier frequency is 
approximately 1575.42 MHz, while the second is approximately 1227.60 MHz. The 
second carrier frequency is predominantly used for military applications. 

Each satellite uses two PRN codes to modulate the first carrier signal. The 
first code is a coarse acquisition (C/A) code. This code is a sequence of 1 023 
elements, called chips, and it modulates at an approximately 1MHz rate. The second 
code is a precise (P) code, which repeats on a seven-day cycle and modulates at a 
10MHz rate. Different PRN codes are assigned to different satellites in order to 
distinguish GPS signals transmitted by different satellites. 
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The navigation data is superimposed on the first carrier signal and the PRN 
codes. The navigation data is transmitted as a sequence of five frames. Each frame is 
six seconds long, and it contains a time stamp for when the frame was transmitted. 
The navigation-data frames also provide information about the satellited clock errors, 
5 the satellite's orbit (z.e., ephemeris data) and other system status data. A GPS satellite 
receives its ephemeris data from monitoring stations that monitor ephemeris errors in 
its altitude, position, and speed. 

Based on the signals transmitted by the GPS satellites, GPS techniques 
typically estimate the location of a GPS receiver by using a triangulation method, 

10 which typically requires the acquisition and tracking of at least four satellite signals at 
the 1.57542 GHz frequency. GPS signal acquisition often involves correlation 
calculations between the received GPS signal and the C/A code of each satellite at 
various phase offsets and Doppler-shifted frequencies. Code phase and phase delay 
are other names for phase offset. For each satellite, GPS acquisition techniques record 

15 the largest-calculated correlation value as well as the phase offset and Doppler-shifted 
frequency resulting in this value. GPS acquisition techniques then identify several 
satellites (e.g., the four satellites) that resulted in the highest-recorded correlation 
values. 

After signal acquisition, a signal tracking process decodes the signals from the 
20 identified satellites at the phase offsets and Doppler values associated with the 
recorded correlation values for these satellites. Specifically, the signal tracking 
process uses the identified phase offset and Doppler values for the identified satellites 
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to extract each identified satellite's navigation data. Part of the extracted data is the 
time stamp information. 

From the extracted time data, the signal tracking process can compute the 
distance between the receiver and the identified satellites. In particular, a satellite's 
signal-transmission delay (i.e. , the time for a signal to travel from the satellite to the 
receiver) can be calculated by subtracting the satellite's transmission time (i.e. , the 
satellite's extracted time stamp) from the time the receiver received the satellite's 
signal. In turn, the distance between the receiver and a satellite can be computed by 
multiplying the satellite's signal-transmission delay by the speed of light. The 
estimated or exact distance between the receiver and a satellite is often referred to as 
the satellite's pseudorange. 

After signal tracking, a triangulation process typically computes the location of 
the GPS receiver based on the computed pseudoranges and the locations of the 
satellites. The location of each satellite identified during signal acquisition can be 
calculated from the satellite's ephemeris data. Theoretically, triangulation requires the 
computation of pseudoranges and locations of only three satellites. However, 
triangulation methods often use the pseudoranges and locations of four satellites 
because of inaccuracies in the receiver clock. 

Some GPS systems also improve their accuracy by using differential GPS 
techniques. Such a technique requires the operation of differential GPS receivers at 
known locations. Unlike regular GPS receivers that use timing signals to calculate 
their positions, the differential GPS receivers use their known locations to calculate 
timing errors due to the signal path. These differential GPS receivers determine what 
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the travel time of the GPS signals should be, and compare them with what they 
actually are. Based on these comparisons, the differential GPS receivers generate 
"error correction" factors, which they relay to nearby GPS receivers. The GPS 
receivers then factor these errors into their calculation of the transmission delay. 

5 Signal tracking has a number of disadvantages. For instance, it is 

computationally intensive, and hence quite time consuming. Also, signal tracking at 
worst needs 6 seconds, and on average needs 3 seconds, of GPS data to extract a 
selected satellite's time stamp, as the time stamp is embedded in the navigation-data 
frames that are 6 seconds long. This requirement, in turn, reduces the speed of the 
10 location-determination process. Finally, signal tracking and reliable decoding of data 
require higher signal power than that needed for acquisition. This becomes a 
significant impediment when the GPS signal is attenuated, as in indoor or urban 
environments. 

Therefore, there is a need for a global positioning system that can quickly 
1 5 identify a GPS receiver's location. There is also a need for a global positioning system 
that only needs a small amount of received GPS data to identify the GPS receiver's 
location. More generally, there is a need for a location-determination system that 
addresses some or all of the above-mentioned needs. 
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SUMMARY OF THE INVENTION 

Some embodiments of the invention provide a location-determination system 
that includes several transmitters and at least one receiver. Each transmitter transmits 
a signal that includes a unique periodically-repeating component, and the receiver 
5 receives a reference signal. Based on the received reference signal, the location- 
determination system identifies an estimated location of the receiver as follows. For 
each transmitter in a set of transmitters, the system computes a phase offset between 
the received reference signal and a replica of the transmitter's periodically-repeating 
component. The system also identifies an approximate location of the receiver and an 
10 approximate receive time for the received signal. The system then uses the identified 
approximate location and time, and the computed phase offsets, to compute 
pseudoranges for the set of transmitters. Finally, the system identifies the estimated 
location of the receiver by using the computed pseudoranges. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The novel features of the invention are set forth in the appended claims. 
However, for purpose of explanation, several embodiments of the invention are set 
forth in the following figures. 

5 Figure 1 illustrates a signal processing circuit that receives a GPS signal and 

generates a number of digital samples from this GPS signal. 

Figure 2 illustrates a process that determines the location of the GPS receiver 
from the digital samples generated by the signal-processing circuit. 

Figure 3 illustrates a location-determination system. 

10 Figure 4 illustrates a process that calculates the approximate location of GPS 

satellites at a particular transmission time. 

Figure 5 illustrates one satellite that is "over" an approximate location and 
another satellite that is not "over" the approximate location. 

Figure 6 illustrates a process that identifies satellites that are over an 
15 approximate location. 

Figure 7 illustrates a process that computes phase offsets for several satellites. 

Figure 8 illustrates the type of code segments that a GPS receiver receives 
from several GPS satellites. 

Figure 9 illustrates a process for computing pseudoranges. 
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Figure 10 illustrates a PRN code received by a GPS receiver, a replica of this 
code stored by the GPS receiver, and a phase-shifted version of this replica. 
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DETAILED DESCRIPTION OF THE INVENTION 

The invention provides location-determination method and apparatus. In the 
following description, numerous details are set forth for purpose of explanation. 
However, one of ordinary skill in the art will realize that the invention may be 
5 practiced without the use of these specific details. For instance, some embodiments of 
the invention are described below by reference to global positioning systems. One of 
ordinary skill will understand that other embodiments of the invention are used in 
other types of location-determination systems. In other instances, well-known 
structures and devices are shown in block diagram form in order not to obscure the 
10 description of the invention with unnecessary detail. 

Some embodiments of the invention are location determination systems that 
determine location of receivers from "reference signals" received by the receivers. As 
used in this document, a reference signal means any type of signal from which 
location information may be derived. Hence, the reference signal can be a GPS 
15 ("global positioning system") signal, a CDMA ("code division multiple access") 

signal, a GSM ("global system for mobile communication") signal, a TDMA ("time 
division multiple access") signal, a CDPD ("cellular digital packet data") signal, or 
any other signal from which location information may be derived. 

In the embodiments described below, however, the reference signal is a GPS 
20 signal that can be used to estimate the location of GPS receivers. On Earth, a GPS 
receiver typically receives a GPS signal that is a composite of several signals 
transmitted by GPS satellites that orbit the Earth. The characteristics of such GPS- 
satellite signals were described in the background section. 
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The GPS receiver can be a standalone device, can be part of another mobile 
device (e.g., a personal digital assistant ("PDA"), wireless telephone, etc.), or can 
communicatively connect to another mobile device (e.g., connect to a Handspring 
Visor PDA through its proprietary Springboard). Several such architectures for the 
5 GPS receiver are described in United States Patent Application, entitled "Method and 
Apparatus for Determining Location Using a Thin-Client Device," filed on December 
4, 2000, and having serial number 09/730,324. The disclosure of this application (i.e., 
United Stated Patent Application, entitled "Method and Apparatus for Determining 
Location Using a Thin-Client Device," filed on December 4, 2000, and having serial 
1 0 number 09/73 0,324) is incorporated herein by reference. 

Some embodiments estimate the location of a GPS receiver by (1) initially 
generating K-digital samples sj, s& from the GPS reference signal received by the 
receiver, and then (2) using the digitized GPS reference data to estimate the location 
of the GPS receiver. The GPS receiver typically performs the digitization operation. 
15 Also, in some embodiments, the GPS receiver digitizes only a portion of the received 
GPS signal. An example of a signal-processing circuit that a GPS receiver can use to 
generate such samples will be described below by reference to Figure 1. 

Some embodiments estimate the location of the GPS receiver from the 
digitized GPS reference data by using a location-determination process that can be 
20 performed either (1) completely by the GPS receiver, (2) completely by another device 
or computer in communication with the GPS receiver, or (3) partially by the GPS 
receiver and partially by another device or computer in communication with the GPS 
receiver. 
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In some embodiments, the location-determination process computes 
pseudoranges for a set of satellites without performing signal tracking on the received 
GPS reference data. In some of these embodiments, the process computes for each 
satellite in the set a phase offset between the samples of the received reference signal 
and a replica of the satellite's PRN code. The system also identifies an approximate 
location of the receiver and an approximate receive time (e.g., an approximate time 
for the start of the received signal). 

The system then uses the receiver's approximate location, approximate receive 
time, and the computed phase offsets to compute pseudoranges for the set of satellites. 
To compute a satellite's pseudorange, some embodiments use the satellite's computed 
phase offset and the receiver's approximate location to calculate the length of all or 
most of the PRN codes transmitted by the satellite from the time that it transmitted the 
signal resulting in the first sample s/ to the time the receiver generated the first sample 
s/. Finally, the system identifies the estimated location of the receiver by using the 
computed pseudoranges. 

In certain circumstances, the estimated receiver location matches the exact 
receiver location. In other circumstances, the estimated receiver location matches the 
exact receiver location to such a high degree of accuracy that it is indistinguishable 
from the exact location to an observer. In yet other situations, however, the estimated 
location differs from the actual location of the GPS receiver by a certain error amount; 
in these situations, some embodiments take steps to ensure that this error (between the 
estimated and actual receiver locations) is tolerable for the particular iocation- 
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determination application. Several more specific location-determination processes 
will be explained by reference to Figures 2-10. 

I. GENERATING THE DIGITAL REFERENCE DATA 

Figure 1 illustrates a signal processing circuit 100 that receives a GPS signal 
and generates K-digital samples from the GPS signal. As shown in this figure, the 
signal processing circuit 100 includes a GPS antenna 105, a GPS tuner 1 10, a clock 
1 15, a down-converter 120, and an analog to digital ("A/D") converter 125. 

The GPS antenna 105 receives a GPS signal s(t) , which on Earth is a 
composite of noise and several signals transmitted by several GPS satellites that orbit 
the Earth. The antenna 105 and its associated circuitry are configured to receive the 
reference GPS signal s( t) at a GPS carrier frequency, which currently is around 1 .57 
gigahertz ("Ghz"). The RF tuner 1 10 receives the GPS signal s(t) from the antenna 
105. This tuner 1 10 is tuned to capture signals at the approximate freqtiency of the 
GPS signal. Hence, the tuner captures the GPS reference signal s( t ) received by the 
antenna 105. 

The RF tuner communicatively couples to the clock 1 15 to receive a clock 
signal. The clock 115 generates one or more clock signals to synchronize the 
operation of the components of the signal-processing circuit 100. This clock also 
receives a synchronizing clock signal 130. The clock 115 can receive the 
synchronizing signal 130 from a variety of sources. For instance, in some 
embodiments, the signal-processing circuit 100 includes an RF processing circuit that 
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(1) captures a radio signal with the synchronizing signal, and (2) supplies this signal to 
the clock 115. 

The synchronizing signal allows the clock to set initially its internal time. It 
also allows the clock to synchronize with a reference clock. In some embodiments, 
the receiver is part of a network (e.g. , a wireless telephone network) and its clock 115 
is synchronized with the network clock so as to keep the receiver clock within a few 
micro- or milli-seconds off the GPS clock. Some embodiments synchronize the clock 
1 15 to be within 1 millisecond or less of the GPS clock. The degree of the receiver 
clock ! s inaccuracy will depend on (1) source and manner of obtaining the 
synchronizing signal 130, (2) the accuracy of the clock of the synchronizing-signal 
source, and (3) the accuracy of the GPS-receiver clock. 

It should be noted that other embodiments use other techniques to maintain the 
receiver clock within a desired time interval of GPS time. For instance, some 
embodiments might have the receiver estimate and compensate for its clock drift 
internally. Other embodiments use high-precision clocks in their receivers. 

As shown in Figure 1, the down converter 120 receives the tuner's output (i.e. 9 
receives the captured GPS reference signal s(t) ). The down converter 120 
transforms the captured GPS reference signal to an intermediate frequency ("IF") 
reference signal s( t) . To do this, the down converter includes in some embodiments 
an IF mixer that converts the frequency of the captured GPS signal to an IF frequency, 
such as 50Mhz. In some embodiments, the down converter also includes one or more 
band pass and amplification stages to filter and amplify the input and/or output of the 
IF mixer. 

12 
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The signal-processing circuit 100 utilizes a down converter so that the A/D 
converter 125 can sample the reference signal at an intermediate frequency as opposed 
to a radio frequency. One of ordinary skill will understand that other embodiments 
can include more than one down converter in their signal-processing circuits. Also, 
5 some embodiments use one or more down converters to convert the GPS reference 
signal to a baseband reference signal, which can then be sampled by the A/D 
converter. 

In some embodiments, the A/D converter's sampling rate is at least twice the 
size of the frequency band, while in other embodiments the sampling rate differs from 
10 this rate. The A/D converter 125 samples the IF reference GPS signal s(t ) that it 
receives from the down converter 120, and outputs K-digital samples s h s K of the 
IF GPS signal s(t) . These samples are from only a portion of the received GPS 
signal, and serve as digital GPS reference data that the location determination 
processes use. 

15 These samples are taken at GPS times tj M ...Jk. The corresponding receiver 

clock times are denoted by ^,... 5 7 K . Hence, 7 X represent the receiver clock time 
recorded when the first sample of the GPS signal is received at the receiver. The first- 
sample time 7 X represents the approximate time for when the first sample was 
generated because it is measured according to the receiver's clock, which is not 

20 completely synchronized with the GPS time in the embodiments described below. 

The A/D converter also outputs the first-sample time 7 X . In this document, a "tilde" is 
generally placed over a variable (like 7 X ) to indicate that the value of the variable is an 

approximation of the actual value of the item represented by the variable. 
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H. ESTIMATING THE LOCATION OF THE GPS RECEIVER 

Figure 2 illustrates a location-determination process 200 that estimates the 
location of the GPS receiver from the digital samples generated by the signal- 
processing circuit 100. In the embodiments described below, a location-determination 
server separate from the GPS receiver performs this process 200. One of ordinary skill 
will also realize that in other embodiments, the process 200 is performed either 
(1) completely by the GPS receiver, or (2) partially by the GPS receiver and partially 
by another device or computer (e.g., a location-determination server) in 
communication with the GPS receiver. 

Figure 3 illustrates a location-determination server 300 that performs the 
process 200. In some embodiments, this server is just one computer, while in other 
embodiments several computers form this server. In some embodiments, these several 
computers may be geographically distributed. This server can be a standalone device 
or it can be part of other devices. 

As shown in Figure 3, the location-determination server 300 receives digital 
GPS reference data sj, sk from client devices 305 through one or more base stations 
310. A client device 305 is either a GPS receiver, contains a GPS receiver, or 
.communicatively connects to a GPS receiver. In addition, each base station detects 
GPS-receiver signal transmissions within a particular region, and relays this 
information to the location-determination server. One of ordinary skill will realize 
that the base station can use a variety of communication architectures and networks to 
relay the signals from the GPS receivers to the location-determination server. In some 
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embodiments, the base station is a cell tower that provides wireless communications 
to the client device 305. 

The location-determination server 300 performs the process 200 for a 
particular GPS receiver whenever it receives the digital GPS reference data si, s K 
5 from the particular GPS receiver. In addition to the receiver's reference GPS data, the 
server 300 uses several other data items to perform process 200, in some embodiments 
of the invention. For instance, some embodiments use an approximate location of the 
GPS receiver to compute pseudoranges for several satellites. In the embodiments 
described below, the approximate location is represented in terms of a three- 
10 dimensional vector 1 , which represents a point in an Earth-centered, Earth-fixed 

coordinate system. Some embodiments specify this approximate location to be within 
±150 Km of the actual GPS receiver location. 

Different embodiments identify the approximate location 1? differently. In the 
embodiments described below, the approximate location is the location of the base 

1 5 station 310 that relays the GPS reference data generated by the receiver. Some of 
these embodiments extract the approximate-location information from the signal 
forwarded by the base station. In other embodiments, the server 300 uses the base- 
station's identification to retrieve the base-station location information from a storage 
structure (such as database 315 of Figure 3) that stores the location of all base 

20 stations. Other embodiments determine the approximate location 7 by measuring 
signal strength between the base station 310 and the client device 305, using time- 
difference of arrival (TOO A), angle of arrival (AO A), and/or sector of arrival (SOA) 
methods at the base stations. Yet other embodiments identify the approximate 

15 
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location as the last CPs-receiver's location that was recorded within a predetermined 
time interval of the current first-sample time. One of ordinary skill will realize that 
other embodiments use numerous other approaches to identify the approximate GPS- 
receiver location. 

5 In some embodiments, the data-processing server also needs to have access to 

differential GPS data, satellite clock-correction data, navigation bits, and ephemeris 
data. The server can receive such data from a variety of sources. For instance, as 
shown in Figure 3, the server can receive such data from one or more reference GPS 
receivers 320 through one or more communication networks 325, such as the Internet. 

10 To perform the process 200 for a particular GPS receiver, the location- 

determination server 300 initially computes (at 205) the approximate location of all 

GPS satellites at the time that these satellites transmitted the signals resulting in the 
first sample sj of the GPS reference data. This computation is further described below 
in Section II. A by reference to Figure 4. 

1 5 The process next identifies (at 21 0) N GPS satellites that are over the GPS 

receiver's approximate location 1 . This identification is further described below in 
Section II.B by reference to Figures 5 and 6. After identifying the N overhead 
satellites for the approximate location 1 , the process 200 computes (at 215) phase 
offsets for P satellites, where P is a number equal to or less than N. This computation 

20 is further described below in Section E.C by reference to Figure 7. 

For each of the P satellites, the process then computes (at 220) a pseudorange 
from its computed code phase, the GPS-receiver f s approximate location, and the 
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satellite's approximate location. As further described below, each satellite's 
approximate location is computed based on the approximate first-sample time 7 X , 
which is based on the receiver clock. The computation of the pseudoranges is further 
described below in Section KD by reference to Figures 8-10. 

The process 200 then uses (at 225) the computed pseudoranges to calculate the 
GPS receiver's estimated location. Some embodiments use triangulation techniques to 
calculate the GPS receiver's estimated location from the computed pseudoranges, as 
described below in Section HE. 

A. Computation of Satellite Locations . 

As mentioned above, the location-determination process 200 initially 
computes the approximate locations of the satellites at the time that these satellites 
transmitted the signals that resulted in the first sample sj. The receiver generated this 
first sample at time 7 X . 

Given ephemeris and differential data, a location y' of a satellite i is for all 

practical purposes a deterministic function of the GPS time t. One set of equations for 
deriving a satellite's location from ephemeris and differential data is provided in Table 
2.3 on page 38 of "Understanding GPS Principles and Applications," by Elliott 
Kaplan, Artech House, 1996. The computed satellite location is typically a three- 
dimensional vector (y i e 5H 5 ) that is defined in an Earth-centered, Earth-fixed 
coordinate system. 
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To compute the location of a satellite i, the process 200 needs to first compute 
the time that this satellite transmitted the signal resulting in the first sample. However, 
the exact transmission time is difficult to compute without performing signal tracking. 
Hence, the process 200 computes an approximate transmission time and location for 
each satellite, instead of computing the precise values of these parameters. 

More specifically, the exact transmission time for each satellite i equals the 
actual first-sample time t\ minus the actual signal-transit delay t, for the satellite. The 
signal-transit delay r* is the time between the satellite's transmission of the signal 
resulting in the first sample Sj and the receiver's generation of the first sample. 
Equation (1) illustrates one manner of computing the exact signal transit delay r/of a 
satellite i. 



In this equation, £ is the receiver's exact position at first-sample time tj, y* (ti-%$ is 
the satellite's precise location at exact transmission time (tj-ri), c is the speed of light, 
r f a is delay due to atmospheric conditions, r r is delay caused by analog processing at 
the receiver. 

The location-determination process 200 can determine atmospheric-delay 
r° by using differential data. It can also retrieve from a storage structure (such as data 

storage 315) the processing-delay r r associated with each GPS receiver. However, it 
is difficult to ascertain (1) the exact first-sample time t ly (2) the receiver's precise 
location £ at exact first-sample time tj 9 or (3) the satellite's precise location y 1 (ti-Tj) 




(1) 
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at precise transmission time (tj-ri). Hence, it is difficult to solve the Equation (1) to 
obtain an exact value for the signal-transit delay t>. 

Therefore, in some embodiments, the location-determination process 200 (at 
205) computes for each satellite / the approximate signal transit delay r } of the 
satellite. Equation (2) below illustrates one manner of computing an approximate 
signal-transit delay r t . This equation assumes that (1) the first-sample time is 7 } , 

(2) the receiver is located at the approximate location 1 , and (3) the satellite's 
location at transmission time equals its location at the approximate first-sample time 
?• 

c 

As mentioned above, some embodiments use the base-station tower location as the 
approximate GPS-receiver location 1 , and retrieve the tower location from a database 
by using the received cell tower identification. 

For each satellite /, this process can thus compute an approximate transmission 
time (t x -tJ based on the approximate first-sample time ^ and the satellite's 
approximated signal-transit delay r l . Accordingly, for each satellite i, the process 
computes the satellite's approximate location y l (7 } - rj at the computed approximate 
transmission time (t x - rj 9 by using the ephemeris and differential data. As mentioned 
above, one . set of equations for deriving a satellite's location given ephemeris and 
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differential data are provided in Table 2.3 on page 38 of "Understanding GPS 
Principles and Applications," by Elliott Kaplan, Artech House, 1996. 

It can be shown that the approximate satellite location that is computed under 
this approach is no more than few meters off from the true satellite location at the time 
5 of transmission. Approximate location errors of less than 150 km will lead to no loss 
in accuracy, while errors in the receive time are less than four meters per millisecond. 
When receiver are part of a network (such as a cellular network), the receiver clocks 
are much more accurate than one millisecond. In the remainder of this document, the 
symbol y* is used as shorthand for the approximate satellite location y l (t x -T t ). 

10 Figure 4 illustrates a process 400 that uses the above-described approach to 

calculate the approximate location of each GPS satellite i at the approximate time (7^ - 
fj that the satellite transmitted the signal resulting in the first generated sample. This 
process is used by process 200 at 205, in some embodiments of the invention. 

As shown in Figure 4, the process 400 initially selects (at 405) a first GPS 
15 satellite / from a list of GPS satellites that it has. The process then determines (at 410) 
the signal-processing delay r r associate with receiver. Some embodiments store in a 
storage structure the processing-delay r r associated with each GPS receiver that can 
be used with the process 200 and the location-determination server 300. Hence, in 
these embodiments, the process 400 retrieves from the storage structure the 
20 processing-delay r r associated with the particular GPS receiver at issue. 

Next, for the selected satellite, the process computes (at 415) the atmospheric 

delay t, a . This atmospheric delay can be derived from the received differential GPS 
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data and the approximate receiver location 1 . One manner for deriving atmospheric 
delay at a location from received differential GPS data is disclosed in Chapter 8 of 
"Understanding GPS Principles and Applications," by Elliott Kaplan, Artech House, 
1996. 

At 420, the process then computes the selected satellite's approximate location 
j/ (t\) at the approximate first-sample time 7 X . In some embodiments, the process 
computes this location at the time 7 X from the ephemeris and differential data by using 
the set of equations provided in Table 2.3 on page 38 of "Understanding GPS 
Principles and Applications," by Elliott Kaplan, Artech House, 1996. 

After computing the approximate location of the satellite at time 7 X , the 
process 400 computes (at 425) the approximate signal-transit delay z i of the selected 
satellite. This process uses the above-described Equation (2) to compute the 
approximate signal-transit delay r. from the approximate satellite location y l (7 X ) 

(computed at 420), the signal-processing and atmospheric delays r r and r* 
(respectively obtained at 410 and 415), and the approximate receiver location 1 . 

Next, the process computes (at 430) the approximate transmission time (7 X -tJ 
for the selected satellite by subtracting the approximate signal-transit delay r t 
(computed at 425) from the approximate first-sample time 7 X . The process then 
computes (at 440) the selected satellite's approximate location y* at the approximate 
transmission time (7 X - r f ) computed at 435. In some embodiments, the process 
computes this location from the ephemeris and differential data by using the set of 
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equations provided in Table 2.3 on page 38 of "Understanding GPS Principles and 
Applications," by Elliott Kaplan, Artech House, 1996. 

After computing the selected satellite's location at the approximate 
transmission time, the process determines (at 440) whether the selected satellite is the 
last GPS satellite on its GPS-satellite list. If not, the process selects (at 445) another 
GPS satellite from this list, and repeats 415 to 435 in order to determine this newly- 
selected satellite's approximate location. Otherwise, the process end as it has 
computed location of all the GPS satellites at the approximate time they transmitted 
their signals that resulted in the first generated sample. 

B. Identification of Overhead Satellites , 

After identifying the approximate location of the GPS satellites, the process 
200 identifies (at 210) the satellites that are currently "overhead." In some 
embodiments, the process 200 identifies the "overhead" satellites by making a 
simplifying assumption that satellites are overhead relative to the client device 305 if 
and only if they are overhead relative to the approximate location £ . 

Hence, in these embodiments, the process 200 examines each satellite j and 
determines whether the satellite is over the approximate location 1 . In some 
embodiments, the process specifies a satellite / as an overhead satellite if and only if 
I'l <!'$*. (Note that the "prime" denotes vector transposition.) In other words, the 
process designates satellite i as an overhead satellite if and only if the inner product of 
the satellite's approximate position vector s t with the approximate location vector 1 is 
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greater than or equal to the inner product of the approximate location vector with 
itself. 

This designation criterion essentially determines whether the magnitude of the 
projection of a satellited approximate position vector y l onto the approximate 
5 location vector i is greater than or equal to the magnitude of the approximate 
location vector £. If so, the satellite is an overhead satellite. 

Figure 5 pictorially illustrates this designation criterion. This figure presents a 
GPS receiver 505 at the approximate location 1 on Earth 510, and two satellites 515 
and 520 that orbit the Earth. The satellite 515 is not an overhead satellite as it is 

10 beneath the horizon 525 of the location 1 . Accordingly, the magnitude of the 

projection of its position vector y 515 onto the approximate location vector t is less 
than the magnitude of the approximate location vector i . On the other hand, the 
satellite 520 is an overhead satellite as it is above the horizon 525 of the location I . 
Hence, the magnitude of the projection of its position vector y 520 onto the approximate 

15 location vector 1 is greater than the magnitude of the approximate location vector i . 

Figure 6 illustrates a process 600 that uses the above-described approach to 
identify the overhead satellites. This process is used by process 200 at 210, in some 
embodiments of the invention. As shown in Figure 6, the process 600 initially 
computes (at 605) the Loner product of the approximate location vector t with itself. 

20 The process 600 next selects (at 610) a GPS satellite from a list of GPS 

satellites that it has. The process then computes (at 61 5) the inner product of the 
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approximate location vector 1 with the selected satellite's approximate location 
vector y 1 , which was computed at 205. 

At 620, the process 600 determines whether the inner product generated at 61 5 
is greater or equal to the inner product generated at 605. If not, the process transitions 
to 630, which will be described below. Otherwise, the process designates (at 625) the 
satellite selected at 610 as an overhead satellite, and then transitions to 630. 

At 630, the process determines whether it has examined all the GPS satellites 
on its list of GPS satellites. If not, the process returns to 610 to select another GPS 
satellite from this list, and repeats the above-described operations in order to 
determine whether the newly-selected satellite is an overhead satellite or not. 

Once the process determines (at 630) that it has examined all the GPS 
satellites, the process ends. Typically, by the time that the process 600 ends, this 
process has identified some integer number N of the GPS satellites as overhead 
satellites. 

C. Computation of Phase offsets . 

After identifying the overhead satellites for the approximate location t , the 
process 200 computes (at 215) phase offsets for the overhead satellite. One method 
for computing phase offsets is described below. However, one of ordinary skill will 
realize that a variety of other methods may be used to compute phase offsets. 

A phase offset for a satellite is the phase difference between (1) the satellites 
PRN code that is part of the received GPS reference data, and (2) a replica of the 
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satellite's PRN code that is used by process 200. The process 200 computes a 
satellite's phase offset by computing an ambiguity function A % (^, f D ) , which is the 
square of the magnitude of the inner product Ri(<f>,f D ) between the received GPS 
reference data with the replica of the satellite's PRN code at various candidate phase 
offsets 0 and Doppler-shift values f D . Equations (3) and (4) provide mathematical 
representations of the ambiguity and correlation functions. 

A,(lfo) = \R,(lf D )\ 2 (3) 
R,(lfo) = ^ - X -0)s{° (4) 

A. k~J K 

In these equations, (1) denotes the GPS time of the Ji h sample, (2) is the fd h 

sample after the candidate Doppler-shift value has been applied to it, (3) the candidate 
Doppler-shift value f D represents the sum of the satellite Doppler and the clock 
Doppler values, and (4) x k (t) is a bandpass filtered version of the satellite's CIA code 
replica x i (t) (modulated by the navigation bits) that the process 200 uses. The 

Doppler-shift value is applied to the A^ 1 sample by multiplying this sample by 
e -2ni<f!F+fD)(tt-h) m 

In some embodiments, the search through candidates ($ 9 f D ) is exhaustive, 
and the phase offset of the candidate pair that resulted in the maximum ambiguity 
value for the satellite is selected as the phase offset of the satellite. As described 
below by reference to Figure 7, some embodiments take this selection a step further. 

25 



WO 02/46786 PCT/USO 1/45828 

These embodiments use an ambiguity threshold ratio to detect useful phase offsets, 
and disregard satellites without useful phase offsets. For instance, some embodiments 
require a satellite's maximum ambiguity to exceed the average ambiguity by a 
constant ratio in order for its phase offset to be considered valid. In other 
5 embodiments, a variety of threshold tests can be used, such as ratio between highest 
ambiguity peak and next highest ambiguity peak. 

Also, rather than computing the inner-product function R t (<j>,f D ) as a sum of 
a series of multiplications for each ^ and f D , some embodiments view i?(-,/ D ) as a 
correlation. In one embodiment, the correlation theorem and fast Fourier transforms 
1 0 ("FFT") are used to improve computation time. 

Figure 7 presents a process 700 for computing useful phase offsets for several 
of the overhead satellites. The process 200 performs (at 215) the process 700 after it 
has identified the overhead satellites. The process 700 initially (at 705) selects an 
overhead satellite and defines this satellite's maximum ambiguity value to be equal 
15 toO. 

Next, the process 700 enumerates (at 710) a number of candidate phase offsets 
0's and Doppler-shift values *s for the selected satellite. The process then selects 
(at 715) a candidate phase offset 0 and Doppler-shift value f D , and computes (at 
720) the ambiguity value for the selected candidate pair by using Equations (3) and (4) 
20 above. 

The process next determines (at 725) whether the ambiguity value computed at 
720 is greater than the maximum ambiguity value recorded for the selected satellite 
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thus far. If not, the process transitions to 735, which will be described below. If so, 
the process (at 730) sets the maximum ambiguity value equal to the value computed at 
720, and stores the satellite's phase offset as the best phase offset thus far. The 
process then transitions to 735. 

5 At 735, the process determines whether it has examined all the candidate pairs. 

If not, the process returns to 715 to select another candidate pair. Otherwise, the 
process determines (at 740) whether the selected satellite's maximum ambiguity value 
exceeds a threshold ambiguity value. If so, the process adds (at 745) the selected 
satellite to the list of satellite's for which a pseudorange needs to be computed, and 
10 records the best phase offset stored at 730 for this satellite as its phase offset. 

From 745, the process transitions to 750. The process also transitions to 750 if 
it determines (at 740) that the selected satellite's maximum ambiguity value does not 
exceed the threshold ambiguity value. At 745, the process determines whether it has 
examined all the overhead satellites. If not, the process returns to 705 to select another 
15 satellite. Otherwise, the process ends. 

D. Computation of Pseudoranges . 

After computing phase offsets for several overhead satellites, the process 200 
then uses (at 220) the receiver's approximate location and the computed phase offsets 
to compute pseudoranges for the set of satellites. To compute a satellite's 
20 pseudorange, the process 200 calculates the length of all or most of the PRN codes 
transmitted by the satellite from the time that it transmitted the signal resulting in the 
first sample si to the time the receiver generated the first sample si. 
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As mentioned above, signal-transit delay refers to this time duration (z.e., it is 
the time between the satellite's transmission of the signal resulting in the first sample 
sj and the receiver's generation of the first sample Sj). During the i* satellite's signal- 
transit delay, satellite i transmits several PRN codes. 

5 Figure 8 illustrates that these PRN codes can be grouped in three categories. 

The first category is the PRN code segment 805 that contains the signal resulting in 
the first sample sj. In the discussion below, this PRN code segment is referred to as 
the first PRN code segment. The first PRN code segment 805 is said to be a complete 
or practically complete segment when the signal resulting in the first sample is part of 
10 the first PRN-code chip, and it is said to be a partial segment when the signal resulting 
in the first sample is part of a chip after the first PRN-code chip. 

The second category includes the full PRN codes 810 that the satellite 
transmitted after the first PRN code but before the first-sample time. In the discussion 
below, these PRN codes are referred to as the full PRN codes. The third category is 

15 the PRN code 815 that the satellite is transmitting when the receiver generates the first 
sample. In the discussion below, this PRN code segment is referred to as the last PRN 
code segment 815. The last PRN code segment 815 is said to be a complete or 
practically complete segment when at the first-sample time the satellite is transmitting 
the last chip of this code, and it is said to be a partial segment when at the first-sample 

20 time the satellite is transmitting a chip before the last chip of this code 

Equation (5) below illustrates one manner of computing the estimated distance 
between the GPS receiver and the f h satellite by summing the length first, last, and 
full PRN code segments after accounting for differential corrections and the satellite- 
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clock's error. In other words, Equation (5) provides the mathematical representation 
of a pseudorange p t for the f h satellite. 



In this equation, (1) p } is the length of the first code segment, (2) m f is the number of 
full PRN codes, (3) / is the full-PRN-code length, which is about 300 Km, (4) v, is 
the i th satellite's transmission phase, (5) c is the speed of light, (6) cv t is the length of 
the last PRN code without adjustment for the satellite clock error, (7) d* is the 
satellite-clock-error correction factor, and (8) d x is the atmospheric-error-correction 
factor. 

_ © 
Accordingly, in Equation (5), (1) the first term (p t ) accounts for the length of 

the first code segment, (2) the second term (mj ) accounts for the length of the full 

PRN codes, (3) the third term ( cv i + d*) accounts for the length of the last PRN code 

after adjustment for satellite-clock error, and (4) the fourth term (d t ) accounts for 

atmospheric delay. 

Figure 9 illustrates a process 900 for computing a satellite's pseudorange 
according to Equation (6) below. 



Equation (6) provides the mathematical representation for computing another 
pseudorange p f for the i° l satellite. Equation (6) is similar to Equation (5) except that 



p l =p l +m s l + (cv i +d t *) + d l 



(5) 



p f = Pi + m f l + d* 4- d t 



(6) 
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it does not include cv, , which is the last-PRN-code's length without adjustment for 
satellite clock error. As further described below, the process 900 does not compute 
cv, , but rather leaves it to a subsequent triangulation operation to account for this 
term for all the P satellites together. 



The length p f of the first code segment 805 of satellite i is the code length 

associated with the phase offset $ of this satellite. Accordingly, this code length is 

referred to as phase length in the discussion below. Figure 10 illustrates the 
relationship between the phase length and phase offset of satellite L This figure 
presents three codes. The first code 1005 is the i th satellite's PRN code in the received 

signal. (In this figure, the illustration of the received PRN code disregards 
degradation and distortion in the received signal.) 

The second code is a replica 1010 of the i th satellite's PRN code that the 
location-determination process uses. Figure 10 illustrates the length of complete code 
replica. The third code is a version 1015 of the replica that has been shifted by a 
phase offset to perfectly match the received code 1005. As shown in Figure 10, the 
phase length associated with the phase offset equals the length of the first code 
segment in the received signal. 

When the phase offset <f> t is represented as a normalized-time value between 0 
to 1, the i th satellite's phase length p t is represented by Equation (7) below: 



1. 



Computing the length of the first PRN code segment 



Pi = fa • 



(7) 
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The phase length can be expressed differently. For instance, it can be expressed 
differently when the phase offset is represented differently. 

2. Computing the length of full PRN code segments 

To compute the length of the full PRN codes 810 that the i th satellite 
transmitted after the first PRN code 805 but before the first-sample time, some 
embodiments first compute the number m l of these complete codes. Equation (8) 
provides a mathematical representation for calculating this number. 



m l = round 



(8) 



As illustrated by this equation, the number of full PRN codes can be estimated based 
on the i th satellite's approximate location y l , the GPS -receiver's approximate location 
£ , the i th satellite's computed phase length p t , and the PRN-code-length /. 

Specifically, this number can be estimated by (1) subtracting the satellite's 
approximate location from the GPS-receiver's approximate location to obtain an 
approximate distance between the satellite and receiver, (2) subtracting the phase 
length plus distance correction factor d* from this distance, (3) dividing the resulting 
value by the PRN-code-length /, and (4) rounding the result of the division to the 
closest integer value. The calculation of the distance correction factor d* is 
discussed in the following subsection. Once the number of full PRN codes is 
estimated, the length of these codes can be obtained by multiplying this number by the 
PRN-code-length /, as mentioned above. 
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3. 



Computing the length of the last PRN code segment and 
adjusting for satellite-clock error. 



So long as the GPS-satellite clocks are all synchronized by accounting for the 
satellite-clock errors, the length of the last code segment 815 is equal for all P 
satellites. Hence, as mentioned above, the process 900 below sets the unadjusted last- 
PRN-code length (cv, ) to zero, and leaves it to a subsequent triangulation operation to 
account for this term for all the P satellites together. 

As mentioned above, Equation (6) includes a distance correction factor d* to 

account for the satellite-clock errors. Equation (9) provides one formula for 
calculating this error correction factor. 



In this equation, s i is the satellite-clock error, which in some embodiments is 

retrieved by the location-determination process from a GPS receiver 320 through the 
communication networks 325. 



As mentioned above, Equation (6) includes a differential correction factor d f . 

This correction factor is to account for signal-path delays due to atmospheric 
conditions and other effects. This differential correction factor can be obtained from a 
variety of sources, including third-party vendors, publicly available data via the 
communication network 325 (e.g. www.ngs.noaa.gov/ CORS/cors-data.html), readings 
of an independent GPS receiver and other sources. 




(9) 



4. 



Adjusting for differential corrections. 
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5. Pseudorange-computing process 900. 

The pseudorange-computing process 900 of Figure 9 will now be described. 
This process initially selects (at 905) one of the P satellites that the process 700 
selected for computing their pseudoranges. Next, the process uses (at 910) the above- 
described Equation (7) to compute the selected satellite's phase length p t from its 
phase offset 0, , which was computed by process 700. 

The process then uses (at 915) Equation (9) to compute distance correction 
factor d* due to the satellite-clock error. As mentioned above, some embodiments 
retrieve the satellite-clock error s g from a GPS receiver 320 through the 
communication networks 325. At 920, the process uses Equation (8) to compute the 
number of the full PRN codes 810 that the selected satellite transmitted after the first 
PRN code 805 but before the first-sample time. The process 900 next computes 
(at 925) the length of these full PRN codes by multiplying the number computed at 
915 by the PRN-code-length /, which is about 300 Km. At 930, the process then 
identifies an atmospheric-delay correction factor d t . 

The process then computes (at 935) the selected satellite's pseudorange by - 
summing the values computed at 910, 920, 925, and 930. The process 900 then 
determines (at 940) whether it has generated pseudoranges for all the P satellites 
identified by process 700. If not, the process returns to 905 to select the next satellite, 
and repeat 910-935 to compute the newly-selected satellite's pseudorange. Otherwise, 
the process ends. 



33 



WO 02/46786 



PCTAJS01/45828 



E. Triangulation of Pseudoranges to Determine Location . 

After computing pseudoranges, the process 200 then uses (at 225) the 
computed pseudoranges to calculate the GPS receiver's estimated location. Some 
embodiments use triangulation techniques to calculate the GPS .receiver's estimated 
location from the computed pseudoranges. 

There are multiple triangulation approaches, such as those discussed on pages 
44-48 of "Understanding GPS Principles and Applications," by Elliot Kaplan, Artech 
House, 1996. One skilled in the art will realize that any one of these approaches can 
readily be used in conjunction with the invention to compute receiver location from 
the computed pseudoranges. 

In some embodiments, the triangulation algorithm uses four satellites to 
compensate for a common uncertainty in the satellite distances, instead of 
compensating for the receiver clock error, which is the traditional GPS interpretation 
for using four satellites during tri angulations. As discussed in Section II.D and Figure 
8, this uncertainty is the last PRN Codes 815. This term is common to all the 
satellites, and thus, by including four satellites, the triangulation algorithm can solve 
for this term. The receiver clock error and the uncertainty of the last PRN Code are 
equivalent manifestations of the same entity. 

One of ordinary skill will realize that the above-described embodiments have 
several advantages. For instance, these embodiments do not perform signal tracking 
to compute the pseudoranges of satellites. Accordingly, these embodiments are faster 
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than some of the prior techniques. In addition, these embodiments do not need several 
seconds of GPS data, as they do not extract the satellite time stamps. 

While the invention has been described with reference to numerous specific 
details, one of ordinary skill in the art will recognize that the invention can be 
embodied in other specific forms without departing from the spirit of the invention. 
For instance, even though the embodiments described above use the commercial GPS 
frequency, other embodiments use the military GPS frequency, while others use both 
the military and commercial frequencies. Thus, one of ordinary skill in the art will 
understand that the invention is not to be limited by the foregoing illustrative details, 
but rather is to be defined by the appended claims. 
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CLAIMS 

We claim: 

1 . A global positioning method comprising: 

a) using a global positioning system ("GPS") receiver to receive a 

5 GPS signal; 

b) computing a pseudorange for a GPS satellite without 
performing signal tracking on the received GPS signal; and 

c) using the computed pseudorange to identify an estimated 
location of the GPS receiver. 

10 2. The global positioning method of claim 1 further comprising: 

a) computing pseudoranges for at least two other GPS satellites 
without performing signal tracking on the received GPS signal; and 

b) using the computed pseudoranges to perform triangulation to 
identify the estimated location of the GPS receiver. 

15 3 . A method of computing a pseudorange for a GPS receiver and a GPS 

satellite, wherein the GPS receiver receives a GPS signal and the GPS satellite 
repeatedly transmits a code, the method comprising 

a) identifying a phase offset between the received GPS signal and 
a replica of the satellited code, wherein the phase offset corresponds to a particular 
20 segment of the code; 
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b) identifying an approximate GPS-receiver location; 

c) identifying an approximate time for the received GPS signal; 

d) computing the pseudorange for the GPS satellite based on the 
computed phase offset, the approximate GPS-receiver location, and the approximate 
time. 

4. The method of claim 3, wherein the approximate GPS-receiver location 
is within 1 50 Km of the actual receiver location. 

5. The method of claim 3, wherein the receiver has a clock that maintains 
time, the method further comprising maintaining the receiver-clock time to be within 
a predetermined time interval of GPS time. 

6. The method of claim 5 5 wherein maintaining the receiver-clock time 
comprises synchronizing the receiver-clock time to match time generated by a 
reference clock. 

7. The method of claim 6, wherein synchronizing the receiver-clock time 
comprises periodically synchronizing the receiver-clock time to match time generated 
by the reference clock. 

8. The method of claim 5, wherein the predetermined time interval is 1 
millisecond. 

9. A method of computing a pseudorange for a receiver and a transmitter, 
wherein the receiver receives a reference signal starting at a particular time and the 
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transmitter transmits a signal with a periodically-repeating component, the method 
comprising 

a) identifying a phase offset between the received reference signal 
and a replica of the transmitter signal's periodically-repeating component, wherein the 
phase offset corresponds to a particular segment of the periodically-repeating 
component; 

b) computing the length of the particular segment; 

c) computing the length of complete periodically-repeating 
components that the particular satellite has transmitted by the particular time after 
transmitting the particular segment; 

d) calculating the pseudorange based on the computed lengths. 

10. The method of claim 9 3 wherein calculating the pseudoranges 
comprises summing the computed lengths. 

1 1 . The method of claim 9 5 wherein identifying the phase offset includes 
correlating the received signal with a replica of the transmitter signal's periodically- 
repeating component at various phase offsets. 

12: The method of claim 9, wherein computing the length of complete 
periodically-repeating components comprises: 

a) calculating the number of complete periodically-repeating 
components that the particular satellite transmitted after the particular segment by the 
particular time; and 
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b) multiplying the calculated number by the length of the complete 
periodically-repeating component. 

1 3 . The method of claim 9, wherein the transmitter is part of a GPS 
satellite and the receiver is a GPS receiver. 

14. The method of claim 1 3 , wherein the satellite includes a clock and the 
clock differs from GPS clock by some error amount, the method further comprising: 

a) identifying the satellite's clock error amount; 

b) generating an error-correction factor; 

c) taking the error-correction factor into account in calculating the 

pseudorange. 

1 5 . The method of claim 1 4, 

a) wherein generating an error-correction factor comprises 
multiplying the error amount by the speed of light; 

b) wherein taking the error-correction factor into account 
comprises adding the error-correction factor to the calculated pseudorange. 

16. The method of claim 13, wherein calculating the pseudorange further 
comprises: 

a) calculating an atmospheric-delay-correction factor for the 

satellite; 
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b) taking the correction factor into account in calculating the 

pseudorange. 

17. The method of claim 16, wherein taking the correction factor into 
account comprises adding the correction factor to the calculated pseudorange. 

18. For a global positioning system ("GPS") that includes a plurality of 
GPS satellites and a GPS-receiver, a method of identifying an estimated location of a 
GPS receiver, wherein each particular GPS satellite transmits a signal that includes a 
periodically-repeating code unique to the particular GPS satellite, wherein the GPS 
receiver receives a signal starting at a particular time, the method comprising 

a) for each particular GPS satellite in a set of the GPS satellites, 

1) determining a phase offset between the received signal 
and a replica of the satellite's code, wherein the phase offset corresponds to a 
particular segment of the particular satellite's code; 

2) computing the length of the code segment; 

3) computing the Iengtii of complete codes that the 
particular satellite transmitted after the code segment by the particular time; and 

4) calculating a pseudorange for the particular satellite 
based on the computed lengths; and 

b) identifying the estimated location of the receiver by using the 
computed pseudoranges. 
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19. The method of claim 18, wherein calculating the pseudorange for each 
particular satellite comprises summing the computed lengths for the particular 
satellite. 

20. The method of claim 1 8, wherein identifying the phase offset for each 
particular satellite includes correlating the received GPS signal with a replica of the 
particular satellited code at various phase offsets. 

21. The method of claim 18, wherein computing the length of complete 
codes for each particular satellite comprises: 



satellite transmitted after the code segment by the particular time; and 

b) multiplying the calculated number of complete codes by the 
length of the satellite's code. 

22. The method of claim 18, wherein each satellite includes a clock and the 
clock differs from GPS clock by some error amount, the method further comprising: 

for each particular satellite, 



a) 



calculating the number of complete codes that the particular 



a) 



identifying the particular satellite's clock error amount; 



b) 



generating an error-correction factor; 



c) 



taking the error-correction factor into account in calculating the 



pseudorange for the particular satellite. 



23 . The method of claim 22, 
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a) wherein generating an error-correction factor comprises 
multiplying the error amount by the speed of light; 

b) wherein taking the error-correction factor into account 
comprises adding the error-correction factor to the calculated pseudorange for the 
particular satellite. 

24. The method of claim 1 8, wherein calculating the pseudorange for each 
particular satellite further comprises: 

a) calculating an atmospheric-delay-correction factor for the 
particular satellite; 

b) taking the correction factor into account in calculating the 

pseudorange for the particular satellite. 

25. The method of claim 24, wherein taking the correction factor into 
account comprises adding the correction factor to the calculated pseudorange for the 
particular satellite. 
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